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Abstract 



Aims. We want to understand the emission mechanism of radio emission from air showers to determine the origin of high-energy 
cosmic rays. Therefore, we study the geometry of the air shower radio emission measured with LOPES and search for systematic 
effects between the direction determined on the radio signal and the direction provided by the particle detector array KASCADE. 
, Methods. We produce 4D radio images on time-scales of nanoseconds using digital beam-forming. Each pixel of the image is 
^-H • calculated for three spatial dimensions and as a function of time. The third spatial dimension is obtained by calculating the beam 
^ ■ focus for a range of curvature radii fitted to the signal wave front. We search this multi-dimensional parameter space for the 
fSJ ' direction of maximum coherence of the air shower radio signal and compare it to the direction provided by KASCADE. 

, Results. The maximum radio emission of air showers is obtained for curvature radii being larger than 3 km. We find that the 

• direction of the emission maximum can change when optimizing the curvature radius. This dependence dominates the statistical 
■ uncertainty for the direction determination with LOPES. Furthermore, we find a tentative increase of the curvature radius to lower 

elevations, where the air showers pass through a larger atmospheric depth. The distribution of the offsets between the directions 
of both experiments is found to decrease linearly with increasing signal-to-noise ratio. Significantly increased offsets and enhanced 
signal strengths are found in events which were modified by strong electric fields in thunderstorm clouds. 

Conclusions. We conclude that the angular resolution of LOPES is sufficient to determine the direction which maximizes the observed 
electric field amplitude. However, the statistical uncertainty of the directions is not determined by the resolution of LOPES, but 
by the uncertainty of the curvature radius. We do not find any systematic deviation between the directions determined from the 
. radio signal and from the detected particles. This result places a strong supportive argument for the use of the radio technique to 
' study the origin of high-energy cosmic rays. 
?H ' 

• Key words, acceleration of particles - elementary particles - radiation mechanisms: non-thermal - instrumentation: detectors - 
■ " " ■ methods: observational - techniques: image processing 
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1. Introduction 

Cosmic rays are particles penetrating the Earth's atmo- 
sphere with energies from GeV up to hundreds of EeVs. 
These particles (mostly protons) have been measured di- 
rectly up to energies of PeV. Most cosmic rays up to PeV 
energies are thought to be accelerated in our Galaxy by 
supernova remnants. However, the origin of the particles 
beyond those energies is still unclear. 

Therefore, the study of high-energy cosmic rays is im- 
portant to find their sources and to understand how they 
were accelerated. However, the flux drops with a steep spec- 
trum from about one particle per year per square meter at 
PeV energies, to about one particle per century per square 
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kilometer at energies of a hundred EeV. Because of this low 
flux and a certain interaction in the Earth's atmosphere, 
cosmic rays with energies beyond about a hundred PeV 
are measured indirectly by a cascade of secondary parti- 
cles. Billions of secondary particles are produced in these 
cascades called extensive air showers and hundreds of thou- 
sands of them reach particle detectors on the ground. 

The relative arrival times of these particles allow to de- 
termine the direction from which the primary particle en- 
tered the Earth's atmosphere. However, the identification of 
the source is hampered by the deflection of the primary par- 
ticle by the magnetic field distribution in our Galaxy. The 
higher the energy, the smaller the deflection of the particle; 
and for very high energies larger than 10^^ eV (10 EeV), the 
proton cyclotron radius becomes larger than the dimension 
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of our Galaxy so that its trajectory possibly points back to 
its source. 

Measurement of these cosmic rays over long pe- 
riods can reveal a clustering of events, which would 
enable the identification of their sources and the 
determination of the particle accelerati on process 
ijAuger Collaboration: J. Abraham et al.ll2007l ). 

Cosmic-ray air showers and th eir origin haye been stud- 
ied with particle detectors since lAuger et al. I (Il989h . The 
position of the shower axis on the ground is determine d 
from the lateral distribution of particles ijMaverl Il992t l , 
which decreases with distance from the shower core; and 
the direction is determined on the particle arrival times 
llMaverlll99.'j). 



Falcke et al. I (|2005f ) and lArdouin et aP l|2005h showed 



that radio emission from cosmic-ray air showers can also 
be used to determine the energy of the primary particle. 
This result revived the detection of air showers in the radio 
regim e, whi ch were first o b served in 1964 by iJellev et"an 
||1965( ) and lAllan fc JonesI (|l966f ). The radio emission in 
particle cascades is produced by electrons and positrons 
being defiected in the Earth's magnetic field. Both emit 
beamed coherent radiati on called geosynchrotron emission 
ijFalcke fc Gorhamll2003f l. Since the electromagnetic emis- 
sion is not attenuated by coUisions, it reaches the ground 
and therefore the longitudinal shower development can be 
studied by observing the emission at several distances from 
the shower axis. The peak voltage of the radio emission de- 
tectable on the ground depends ap proximately linearly o n 
the energy of the primary particle ijHuege fc Falckell2005f l. 
In addition, the radio emission increases with increasing ge- 
omagnetic angle betwe en the cosmic-ray t rajectory and the 
Earth's magnetic field (jPalcke et al.ll2005l l. 

A new generation of radio telescopes like LOFAR and 
LOPES can digitize the signals from cosmic-ray air showers 
detected by many simple dipole antennae providing them 
with high time, frequency, and spatial resolution. These 
powerful interferometers can be used to determine the di- 
rection of origin of the beamed and coherent emission pro- 
duced in air showers. 

In this article, we present the determination of the direc- 
tion of several hundreds of cosmic-ray air showers measured 
in the ra dio regime with the LOFAR Prototype Station 
(LOPES, iHorneffer et"ani2004| ) at energies above 10^^ eV. 
We determine the exact direction of maximum emission by 
searching a small 4D radio map in the direction provided 
by the triggering KArlsru he Shower Core and Array parti- 
cle DEtector (KASCADE. lAntoni et £01120031 ). The result- 
ing distance between the radio-direction and the particle- 
direction is analyzed for systematic effects. 

The data acquisition and event selection is described in 
Sect. [21 the direction determination is explained in Sect. 
[3l The dependencies of the shower emission maximum are 
elaborated in Sect. 13. 3[ and the results are discussed in Sect. 

a 



80 Megasamples per second, and written to disk. The dig- 
itally stored time-series for each antenna is processed off- 
line, where we perform absolute calibration, we reduce radio 
frequency interference (RFI) by down- weighting of narrow- 
band lines in frequency, and form a beam by adding up 
or correlating the time-shifted antennae for the direction 
given by KASCADE. In addition to the direction informa- 
tion, KASCADE provides parameters such as the shower 
core position, electron number and muon number allowing 
the estimation of the primary energy. More details on the 
L OPES data, acquisition and data processing can be found 
in lNigl et all (|2007. in prep.t ). 

The shower direction is reconstructed by KASCADE 
with an accuracy of about 0.1 degrees at the zenith and 
the shower core positio n with an accuracy of about 1 me- 
ter ijAntoni et al 1 I2003D . both for an event with an electron 
number of 10^'^ (~ 10^^'^ eV). The directional accuracy im- 
proves with increasing shower size and reduces with inclina- 
tion angle. Beyond 42 degrees zenith angle, the directional 
uncertainty is strongly increased due to the reduced accep- 
tance of the particle detectors. The mentioned accuracy of 
KASCADE is valid for showers hitting the particle detector 
array within a centered circle with a radius of 91 meters, 
so that the core is at least 10 meters inside the array. 

For this work, 664 LOPES events were pre-selected with 
a lower limit on the muon number of 10^'^ (> 10^^ eV) and 
a cut on the distance of the shower core from the array 
centre of 91 meters. These events have been acquired from 
2005-11-16 to 2006-07-23. 



3. Direction Determination 

The phasing of an array by digital beam-forming is usually 
done in direction of a source at infinite distance. However, 
the air shower radio emission is produced in the Earth's 
troposphere by an extended emission-region and therefore, 
the wavefront is not fiat. As a first approximation, a spher- 
ical wavefront was assumed for the radio emission reaching 
the antennae of LOPES. The spherical wavefront is defined 
by a specific curvature radius (beam focus) and a shower 
core position on the ground. The resulting distance of the 
beam focus does not necessarily coincide with the distance 
of the maximum emission in the air shower. This needs to 
be investigated in shower simulations and dedicated mea- 
surements in order to study the real shape of the shower 
front and its relation to the position of the maximum radio 
emission and the position of the maximum number of parti- 
cles, respectively. The relation between the maximum in the 
radio shower and in the particle shower might enable the 
determination of the particle species from the shower radio 
emission only. This topic is not part of this work, however 
Monte Carlo simulations bv iHuege fc Falckd l|2005l . Fig. 5) 
show an indication of a more complex shower wavefront 
than a spherical one. 



2. Data acquisition 

LOPES consists of 30 inverted- V-shaped dipoles placed 
between the particle detector stations of the KASCADE 
experiment. LOPES receives a coincidence trigger from 
KASCADE and stores 0.4 ms of radio data before and af- 
ter its reception. The signal of each antenna is amplified, 
filtered in the band from 40 MHz to 80 MHz, digitized with 



3.1. Four-dimensional radio image 

For each cosmic-ray event, all pixels of a four-dimensional 
radio image were calculated spanning azimuth, elevation, 
curvature radius, and time. The image was calculated with 
an extension of 50 x 50 pixels, with an angular resolution 
of 0.2°. The map was centered on the direction given by 
KASCADE. 
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The shower core position for the spherical beam-forming 
was taken from KASCADE. The curvature radius for the 
beam-forming algorithm was chosen to range from 1 km 
to 10 km in 37 steps of 250 m. The time resolution was 
the LOPES sample-time bin of 12.5 ns. The time window 
was chosen to be 1.6 /xs (128 samples), so much larger than 
the length of a typical cosmic-ray radio pulse of a few sam- 
ples in time, to avoid distortion of the pulse in the image 
due to edge-effects in the beam-forming process. The signal 
has been integrated over the LOPES band from 43 MHz 
to 74 MHz. The shower core position, which is used for 
the beam-forming process remains constant at the value 
provided by KASCADE. The intensity of the cosmic-ray 
event in the image corresponds to the power measured in 
the LOPES beam. The imager performs beam-forming by 
adding up the antennae and not by correlating them, as it 
is used for calculating the E-field amplitude at the position 
of the maximum from the image. 

Oscillations in the signal shorter than the inverse of 
the bandwidth are not real but they are caused by over- 
sampling a bandpass broadened signal. Therefore, the final 
4D image was Hanning-smoothed in time by a weighted 
average of every pixel with its two neighbouring pixels. 

3.2. Emission maximum 

The emission maximum in the 4D image was located, 
firstly, by determining the time of maximum emission 
for each slice in curvature radius by selecting the 2D- 
azimuth/elevation-slice containing the pixel with the max- 
imum intensity. For this search in time, 41 of the 128 
available samples around the window center were scanned. 
Secondly, each of the 2D curvature-slices at the tempo- 
ral maximum is searched for the position of the bright- 
est point sourc e (FIND POINTSOURCES function of 
CASA/AIPS+ +, |wwwTaips2. nrao.edu, vl.9, build 1360, 
IShannorJll996t l. A plot of an example event can be found 
in FiglH In case multiple sources were found, the direction 
of the strongest one was saved and the other sources were 
classified as side- lobes. As a result, the position in time 
and the direction in azimuth-elevation were calculated as a 
function of curvature radius. 

The two-dimensional spatial shape of the cosmic-ray 
radio emission in the image is mainly determined by the 
beam-shape (point spread function, PSF) and depends on 
the layout of the selected antennae. The PSF of the emis- 
sion is assumed to be spherical at the zenith. At lower el- 
evations, the resolution in elevation of LOPES is reduced 
by the sine of the elevation angle. A theoretical uncertainty 
for the deter r ninati on of the direction is obtained following 
iTavlor et all (|1999[ ): 

Aa™„ = r, ^ o (1) 

2 ep D cosO 

Here Aamin is the minimum statistical uncertainty of 
the determined direction, ep is the electric field strength 
of the detected pulse, is the noise in the LOPES beam, 
A is the wavelength, D is the diameter of the LOPES an- 
tenna array, and 9 is the zenith angle. The angular reso- 
lution d = X/Dco&e or beam width of LOPES is 2.2° for 
the zenith (A = 5 m, D cosd ~ 130 m). An example-event 
of LOPES at an estimated primary energy of 50 PeV, an 
electron number of 10^ '', and a zenith angle of 13.0° (for 



LOPES and KASCADE) was obtained with a E-field peak 
signal-to-noise ratio (SNR) of ep/cAr = 8.5 resulting in an 
uncertainty of 0.13°, which is about the same accuracy as 
mentioned for KASCADE. 



3.3. Dependencies of the measured electric field 

In the time domain, the air shower pulse is expected at ap- 
proximately 1.8 /IS be fore the trigger from K ASCADE ar- 
rives at the antennae (|Horneffer et al.l[2004( ). Fig. [T] shows 
the intensity maximum of the image-slices as a function of 
time (x-axis) and curvature radius (different curves). The 
curves for the different curvature radii from 1 km to 10 km 
show that the intensity in the signal at —1.8 /xs is increas- 
ing until the best fit is found, whereas the intensity of the 
background noise stays the same. 

The intensity of the maximum source in each skymap 
increases steeply, when increasing the curvature radii from 
1 km to about 3 km, which is due to increasing coherence of 
the cosmic-ray air shower pulses in the individual antenna 
signals. For larger curvature radii, the intensity either keeps 
increasing asymptotically or a maximum is found, as in Fig. 
m For the plotted event the intensity peaks at 3 km. 

The intensity as a function of azimuth, elevation, and 
curvature radius (different curves) is shown in Fig. [3l The 
different curves show with increasing curvature radius the 
increase in the intensity of the pulse, while the intensity of 
the side-lobes next to the pulse stay the same or are reduced 
in intensity. 

In addition to the change in the intensity, in the plot 
with intensity versus elevation (bottom panel of Fig. [3]), 
the position of the pulse-maximum shifts with increasing 
curvature radius to higher elevations by about two degrees. 
We attribute this to the offset between the phase center of 
the LOPES array and the shower core position as provided 
by KASCADE. The larger the distance between both, the 
more the position of the emission maximum changes, when 
varying the curvature radius for the beam-forming. 

This beam-forming effect was tested by artificially in- 
troducing an extra delay by half a sample to all antenna 
data. This subsample shift is equivalent to a half-sample 
interpolation of the data, when correcting the time-axis by 
half a sample in the opposite shifting direction. The result 
of this method is displayed in Fig. IH which is showing the 
skymaps for the unmodified and for the shifted data, respec- 
tively. It turns out, that the interpolation of the antenna 
data, from a sample time of 12.5 ns to 6.25 ns, can further 
optimize the electric field amplitude of the air shower sig- 
nal found in the beam. This is due to the sampling of the 
LOPES antennae signals in the second Nyquist-zone from 
40 MHz to 80 MHz with a sample frequency of 80 MHz. 
Second Nyquist-sampling preserves all signal information, 
which can be recovered by digital upsampHng with, in our 
case, a factor of two. 

We effectively did this by analyzing the interpolated 
data. In addition to an increase in the electric field ampli- 
tude, the pulse shifts in position, as observed when optimiz- 
ing the curvature radius. Therefore, we conclude that the 
upsampling is necessary for the determination of the correct 
electric field amplitude. As expected, this will increase the 
accuracy of LOPES in field strength and direction estimate. 
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Figure 1. Maximum electric field amplitude of the image 
as a function of time for difiFerent curvature radii. The am- 
plitude of the pulse at —l.Sfis increases with beam-forming 
curvature radii (as shown for another event in Fig. [2|) . 
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Figure 2. Maximum electric field amplitude of the image 
as a function of beam-forming radius (distance from the 
observer to the beam focus) for an event with an estimated 
primary energy of 3 x 10^^ eV. 



4. Event selection 

The 664 pre-selected events (see Sect. ^ were processed 
as described above, including the half-sample shift. Out of 
these 664 events, 232 were selected having a minimum E- 
field peak SNR of 6. The cut on the SNR of the events has 
been appHed since for low SNR the measured noise from 
the KASCADE particle detectors dominates the air shower 
emission, which leads to a wrong fit when searching for the 
pulse profile in the time window of the 4D image. 

The plot in Fig.[T]shows this window and the air shower 
pulse lies at —1.8 ^s, with increased curvature radius rising 
above the particle detector noise of KASCADE at —1.75 /is 
to —1.6 iJ,s. For events with a SNR below 6 the cosmic-ray 
radio pulse does not exceed the emission from KASCADE 
and the wrong signal is fitted. A third selection criterion 
on the coherence of the signal, obtained from the cross- 
correlation beam, and a fourth criterion on the elevation 
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Figure 3. Electric field amplitude as a function of azimuth 
(top) and elevation (bottom) at the time of maximum emis- 
sion in the image for different curvature radii. The ampli- 
tude of the curves increases with beam-forming radii (as 
shown for another event in Fig. [21). The radius of curvature 
is increased by steps of 250 m. 



angle not to be smaller than 15° reduced the number of the 
selected events further to 62. 

The cross-correlation beam is formed by correlating all 
paired combinations of the shifted antenna signals. The 
third selection criterion required the ratio of the electric 
field amplitudes of the cross-correlation beam and the total 
power beam to be larger than 0.73, which is one sigma above 
the average fraction of all events. This selection favors in- 
clined events, since RFI from the KASCADE detectors is 
reduced, due to thei r reduced acceptanc e to larger zenith 
angles (see Fig.lSland lPetrovic et al.ll2Q07h . This way, events 
which consist of RFI but have a SNR above 6 are rejected. 
Such events were found to have small curvature radii, since 
most RFI is emitted locally. We note, that by this criterion 
also events were rejected from air showers which did not 
reach the required level of signal-coherence. 

The rejection of RFI events can be tested by plotting a 
histogram of the distribution of the curvature radii. It can 
be seen in Fig. [6] that indeed all events which had an emis- 
sion maximum at a curvature radius below 2000 m are not 
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part of the filtered distribution, these events are considered 
as RFI events. 

Among the 62 selected events, 17 events have an electric 
field strength of more than twice the expected field strength 
estimated from the muon number, geomag netic angle, and 
shower axis distance (|Horneffer et al.l[2007^ . The measured 
electric field strength of the 62 selected events is plotted 
versus the directional offset in Fig.[7l In this plot, the 17 en- 
hanced events are marked with squares and the 5 encircled 
events were confirmed to be recorded during thunderstorm 
activity with lightnings in the vicinity of LOPES. The error 
bars in the plot give the minimum uncertainty calculated 
with Eq. m 



5. Results and discussion 

The 5 encircled events of Fig. [7] are candidates for en- 
hanced showe r radio emission by electric fields in thunder- 
storm clouds ( Buitink et al.|[2007h . These five thunderstorm 
events show large deviations in position, which suggests 
that the net force of the electric field in thunderstorm clouds 
can change the direction in which the emission is beamed. 
One of the 5 encircled events did not exceed the expected 
field strength by a factor of two, but had a strong offset 
in direction. The other 13 enhanced events could not be 
associated with thunderstorm activity, however, an effect 
by an increased geoelectric field in a cumulonimbus cloud 
without lightning is not excluded. Therefore, the enhanced 
events and the thunderstorm events were not taken into 
account for the following interpretation of the results. The 
remaining 44 events show an approximately linear decrease 
in directional offset with increasing signal strength. 



5.1. Curvature radius 

As mentioned in Sect. 13.31 the intensity detected in the 
pulses increases quickly with curvature radius. In fact, 90% 
of the maximum intensity is obtained for the selected 44 
events at (2600±800) m and the curvature radii at the 
maxima are (7000±2000) m (one sigma statistical uncer- 
tainties). The statistical uncertainty for the radii at the 
maxima is large, because the change in the intensity at 
large curvature radii is small and cannot be determined 
more accurately with an array of the size of LOPES. An 
approximation for the statistical uncertainty of the curva- 
ture radius can be obtained adapting Eq. [TJ 



P J, 



± 



A 



D ep Dcos9 



(2) 



Here Ap is the uncertainty of the radius of curvature 
p = 7000 m and for D = 130 m, A = 5 m, eAr/ep = 1/16, 
and = 40° results in ±1400 m, which is consistent with 
the above measured ±2000 m. 

The histogram of the curvature radii in Fig. [6] shows 
that a few events accumulate at 10 km indicating that 
those events could be fitted with even larger curvature radii. 
However, on average, the intensity does not change signif- 
icantly increasing the curvature radius over the last 2 km, 
and therefore the wavefront is considered to be nearly pla- 
nar. A tendency to larger curvature radii is observed for 
events coming from lower elevation angles, which travelled 
a longer path through the Earth's atmosphere (see Fig. [8]). 



Furthermore, we showed that the variation of the cur- 
vature radius can change the direction of the pulse in eleva- 
tion. We found that the angular offset in elevation from the 
optimum position as a function of curvature radius is on 
average zero for the 44 selected events and therefore does 
not introduce a systematic effect. However, the spread must 
be considered as a statistical uncertainty of ±0.3° for the 
directions determined with LOPES. 



5.2. Upsampling 

Additionally, we showed in Sect. 13. 3^ that shifts in the di- 
rections occur when maximizing the intensity of the pulse 
by interpolating the antenna signals. The change in the 
positions introduced by the half-sample shifts amounts to 
(0.7 ± 0.9)°, for the 44 selected events. This result de- 
termines the improvement of the directional accuracy of 
LOPES. Effectively, the digitization step error was reduced, 
after the signal was upsampled. 

The optimization of the beam-forming by the upsam- 
pling (half-sample shift) did not significantly increase the 
electric field amplitude of the selected events with (1 ± 
3) m~^ MHz~^ compared to th e values determine d 
with the standard analysis software ijBahren et al.l 120061 ). 
However, it reduced the number of events with more than 
one source in the optimized skymap by 30%. Thus, the in- 
tensity of the environmental background noise in the side- 
lobes of the event beams was reduced by the optimized 
beam-forming. 

5.3. Systematics 

Fig. [5] shows arrows pointing from the position recon- 
structed by LOPES to the position reconstructed by 
KASCADE. The arrows for all events are isotropically dis- 
tributed in the circular area that covers the sky. However, 
the arrows of the 44 events (selected on their E-field peak 
SNR) show a clear preference for coming from the North. 
Due to the inclination of the magnetic field of 65° at the 
latitude of LOPES of 49°, the geomagnetic angles are larger 
to the North than to the South and therefore a larger signal 
is produced in the air shower by the geomagnetic effect. 

The components of the arrows in the East- West direc- 
tion and the North-South direction are plotted in Fig. [9l 
Both data sets follow a normal distribution with a full 
width at haff maximum (FWHM) of 1.2° (x^ed = 2-37) 
for the North-South direction and of 0.6° (Xred = 0-88) 
for the East-West direction. These values are consistent 
with the absolute offset of (1.3 ± 0.8)° (see next section). 
Furthermore, both distributions have an average consistent 
with zero and therefore, no systematic offset is found in the 
directions determined by LOPES and KASCADE. 

5.4. Statistics 

The obtained offsets between the directions provided by 
KASCADE and the directions determined from the LOPES 
radio data for the 44 selected events result in (1.3 ± 0.8)°. 
The histogram of these offsets and of all 664 events peak 
at the same bin of 1.5° suggesting that this is the sum of 
the average uncertainty for the direction determination of 
both experiments for the whole range of elevation angles. 
All events having larger offsets are dominated by increased 
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noise from the particle detectors and they are not part of 
the selected 44 events (see Fig[TO|). 

The clear peaks in the curves for all 664 events of Fig. [9] 
and Fig. [To] indicate that the direction was successfully de- 
termined within 1.3° for more than the 44 events. However, 
those events could also consist of RFI. The distinction be- 
tween cosmic-ray events and RFI events could be improved 
by a better knowledge of the position and shape of the air 
shower pulse in time. This knowledge would allow a fit of 
the pulse in Fig. [U which would enable the direction iden- 
tification of events detected with less intensity than the 
corresponding KASCADE noise in the LOPES beam. 

The 44 selected events have an average zenith an- 
gle of 40° and the average electron number detected by 
KASCADE is 10^ ®. The accuracy of KASCADE for an air 
shower with this electron number is ±0.2°. This value is 
corrected for the reduced acceptance of the particle detec- 
tors by dividing with the sine of the zenith angle of 40°. 
Therefore, we conclude that the measured offsets are dom- 
inated by the accuracy of the direction determination with 
LOPES of about ±1.1°. 

This accuracy of LOPES lies above the statistical un- 
certainty determined with Eq. [T]of ±0.2°, since the fit of a 
curved wavefront with a curvature radius p introduces an 
additional uncertainty: 



1 e 



Atttot = ±- 



N 



X 



'2ep D cos9 



(3) 



6. Conclusions 



The 4D mapping of radio emission is a powerful new tool for 
the study of radio sources within the Earth's atmosphere, 
in particular cosmic-ray air showers, because of the infor- 
mation in the third spatial dimension provided by the dis- 
tance of the beam-forming focus. We used an array with 
only 30 single dipoles and succeeded to map out the radio 
emission from air showers, to search for the position on the 
sky, and to fit the curvature of the detected wave front of 
the shower pulse. For this work, we used KASCADE as a 
reference. The mapping-tool can also be used for astrophys- 
ical imaging of fast pulses, since images can be calculated 
on sample time resolution. Furthermore, the third spatial 
dimension of the images allows to distinguish between the 
location of all kinds of transient events occurring inside and 
outside the Earth's atmosphere. 

We found that the position of the emission maximum in 
the map, calculated in the beam-forming process, is sensi- 
tive to the array layout and the phase errors of the antenna 
electronics and therefore changes when beam-forming pa- 
rameters such as the curvature radius or the filtered band 
are varied. We optimized the beam-forming by upsampling 
of the antenna signals, which strongly improved the accu- 
racy of the direction determination. 

We conclude that the angular resolution of LOPES is 
sufficient to localize the air shower axis and to maximize 
the received electric field amplitude. However, the accuracy 
for the direction determination is determined by the uncer- 
tainty of the curvature radius. An exact determination of 
the shape of the shower wavefront seems to be necessary to 
further improve the localization of the shower radio maxi- 
mum in the beam-forming process. Furthermore, the direc- 
tional accuracy could possibly be improved by larger base- 
lines and by calculating every pixel of the searched skymap 



with the cross-correlation beam (sum of all antenna pair 
correlations) and by further suppressing the KASCADE 
noise. Also, more phase-stability of the electronics and a 
larger antenna layout optimized for beam-forming in all 
azimuthal directions is desirable. 

We found that the best fit in the curvature radius in- 
creases for events coming from the horizon after passing 
through a larger atmospheric depth. Furthermore, we found 
an approximately linear decrease in the directional offsets 
between the LOPES positions and the KASCADE positions 
with increasing signal strength. 

Furthermore, we found a few events with field strengths 
enhanced by more than twice the expected field estimated 
from the muon number, geomagnetic angle, and shower axis 
distance. Some of these events were associated with thun- 
derstorms in the vicinity of the telescope, which can en- 
hance the emission through strong electric fields in clouds. 
These events showed the largest offsets between the position 
measured with the radio antennae of LOPES and the par- 
ticle detectors of KASCADE. This very important result 
suggests that the charged particles of air showers experi- 
ence additional deflection by electric flelds in thunderstorm 
clouds. Thus, we measured for the first time a geoelectric 
effect in addition to the geomagnetic effect. 

The offsets between the directions measured with both 
instruments were analyzed for dependencies on the mea- 
sured electric field strength, the curvature radius, the az- 
imuth angle, the elevation angle and the geomagnetic field 
angle. Due to the geomagnetic effect and the geoelectric ef- 
fect discussed here, one could expect a net directional offset 
between the muons, the electrons and the electromagnetic 
emission at the shower maximum. However, the analyzed 
sample of LOPES events did not show any significant sys- 
tematic effect for the mentioned dependencies. 

This result adds another strong argument for the usage 
of the radio detection technique for the study of the arrival 
directions of high-energy cosmic rays. An interesting next 
step will be to investigate with simulations how the dis- 
tance of the beam focus (curvature radius) relates to the 
maximum of the particle shower. 
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Figure 4. Example of radio skyinaps for one LOPES event 
with original sampling (top) and with a half-sample shift 
applied to all antennae signals before beam-forming (bot- 
tom). The split maximum in the top plot is not real and 
merges to one maximum with higher temporal resolution. 
The skymaps show the ai r shower radio emission in s tereo- 
graphic projection (STG. ICalabretta &: Greisenll2002l ) with 
a black dot at the maximum. The second black dot in the 
center of the map and offset from the maximum indicates 
the direction provided by KASCADE. 
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Figure 5. This plot shows an arrow for each event connect- 
ing the position of LOPES and KASCADE projected onto 
the sky with the zenith in the origin. The azimuth angle 
rotates clock-wise starting from North pointing upwards. 
The thin arrows indicate all 664 pre-selected events and 
the thick encircled arrows highhght the 44 selected events. 
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Figure 6. Histogram of the number of events as a function 
of curvature radius. The upper (dotted-blue) distribution 
shows all 664 pre-selected events and the lower (soHd-red) 
distribution shows the selected 44 events. The bin size is 
250 m in curvature radius. 
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Figure 7. Absolute angular separation between the LOPES 
and the KASCADE position as a function of electric field 
amplitude. The encircled events are possibly enhanced in 
field strength in thunderstorm clouds at the time of the ob- 
servation. The events marked with squares have a measured 
field strength of more than twice the expected field strength 
estimated from the r auon number, geomag netic angle, and 
shower axis distance (jHorneffer et al ]l2n07l f. The lower limit 
on the statistical uncertainty in the direction of each event 
is plotted as an error bar (according to Eq. [1]) . 
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Figure 8. Curvature radius as a function of the inverse co- 
sine of the zenith angle for all 62 events. The enhanced and 
thunderstorm events are encircled in this plot. The dashed 
line at 10 km indicates the maximum curvature radius that 
was fitted to the wavefronts. 
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Figure 9. Histograms of the angular offset in the air shower 
core position on the sky at the distance of the curvature ra- 
dius of the air shower emission maximum, as found in the 
LOPES radio signal and as determined with the direction 
provided by KASCADE. The distribution in East- West off- 
set is plotted at the top and in North-South direction at the 
bottom. Each plot shows the distribution for all 664 events 
(dotted) and for the 44 selected events (solid), both fitted 
with a Gaussian (dash-dotted and dashed). 
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Figure 10. Distribution of events as a function of direc- 
tional offset between LOPES and KASCADE. The upper 
(dotted) distribution shows all 664 pre-selected events and 
the lower (solid) distribution shows the selected 44 events 
with the maximum at (1.3 ± 0.8)°. The bin size is 0.5°. 



